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The angular and polarization properties of the photon pair emitted in the two-photon decay of 
2s hyperfine states in hydrogenlike ions are investigated within the relativistic Dirac framework and 
second order perturbation theory. The studied transitions are 2si/2 {F = 1, 0) — >■ lsi/2 {F = 1, 0) in 
Hydrogen atom and 2si/2 {F = 4, 3) —s- lsi/2 {F ~ 4, 3) in hydrogenlike Uranium ion. Two different 
emission patterns are found: For non-spin-flip transitions, the angular correlation (i.e. the angular 
distribution of the emitted photons) is of the type ~ 1 -I- cos^ 6, while the linear polarizations of 
the emitted photons are approximately parallel one another; For spin-flip transitions, the angular 
correlation is of the type ~ 1 — 1/3 cos^ 6, while the linear polarizations of the emitted photons are 
approximately orthogonal one another. 

PACS numbers: 31.10.-l-z, 32.10.Fn 
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I. INTRODUCTION 



The theoretical formalism of the two-photon decay in 
atoms and ions has been introduced by Max Born's PhD 
student Gocppert-Mayer in 1931 [l|. Many aspects of 
such a process, hke the total decay rate and the spec- 
tral distribution, have been extensively investigated in 
the context of few-electron atoms and ions, both in the- 
ory and experiments [a-ll3l|. Recently, some interest has 
been also devoted to the angular and polarization corre- 
lations of the two emitted photons [l4|-[l8ll , following the 
pioneer semi-relativistic work of Au [l9| . Since the be- 
ginning of these studies, two-photon transitions in atoms 
and ions have been not only a challenging process to the- 
oretically and experimentally investigate but also a use- 
ful tool with which to explore different physical areas. 
Already in 1940, for instance, Breit and Teller derived 
that the double photon emission was the principal cause 
of the decay of interstellar hydrogen atoms from their 
metastable 2s state [201, while, more recently, polariza- 
tion correlations of the emitted photons have been em- 
ployed to successfully test violations of the Bell inequality 
[2l| - |24| . Furthermore, two-photon transitions have been 
proposed as a tool to measure weak interaction proper- 
ties [2l,[2i|. 

In the present contribution, the angular and polarization 
correlations of the photons emitted in two-photon decays 
of 2s hyperfine states in hydrogenlike ions are presented. 
By comparing with previous theoretical results, where 
the two-photon decays of bare electron states in hydro- 
genlike ions have been analyzed, our calculations show 
the impact of the nuclear angular momentum (spin) on 
the angular and polarization properties of the emitted 
photons. 
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Nuclear features which are not easily accessible via bare 
nuclear physics might be successfully explored by accu- 
rately studying nuclear effects on atomic processes. A 
brilliant review of such effects, at the borderline between 
atomic and nuclear physics, has been recently written 
by Palffy [231 ■ Along this line, owing to the significant 
progress that has been recently made in the development 
of position-sensitive photon detectors J28l - l31| , the present 
paper may also outline an alternative route to probe nu- 
clear spin properties by using two-photon transitions. 
The paper is structured as follows. In Sect.|n]the theory 
needed for the analysis will be explained. Particularly, in 
Sect. Ill Al the states of the overall -electron plus nucleus- 
system will be defined. In Sect. IIIBl a detailed formu- 
lation of the two-photon bound-bound transition ampli- 
tude will be presented, within the framework of second 
order perturbation theory. Then, in Sect. Ill Cj the differ- 
ential decay rate and the angular correlation of the pho- 
tons will be defined in terms of the transition amplitude. 
In Sect, mil results for the angular correlation function, 
for the transitions 2si/2 (-F = 1,0) — > lsi/2 (i^ = 1,0) 
in Hydrogen atom as well as for 2si/2 {F = 4, 3) ^• 
1^1/2 {F — 4, 3) in hydrogenlike Uranium (gl^U) ion, will 
be presented and discussed. The polarization properties 
of the emitted radiation will be also analyzed in detail 
for a few of these transitions. Finally, a brief summary 
will be given in Sect. IIVI 



II. THEORY 

A. Construction of the overall set of states 

The presence of the nuclear spin has a twofold effect 
on the states of the overall hydrogenlike system. First, 
the energies of the ionic metastable states are slightly 
shifted, mainly due to the magnetic dipole interaction 
that nucleus and electron experience. This energy cor- 



rection can be described by using first order perturba- 
tion theory with additional contributions, such as the 
relativistic, Bohr-Weisskopf, Breit-Rosenthal and QED 
contributions [32l . |33| . Since this energy correction is too 
small to affect considerably the angular and polarization 
properties of the emitted radiation, it will be totally ne- 
glected in the following. 

Secondly, the ionic states acquire a new quantum num- 
ber, usually denoted by F, that represents the total an- 
gular momentum of the overall -nucleus plus electron- 
system. The overall ionic state can be thus described, to 
a first approximation, by coupling the nucleus and elec- 
tron angular momenta: 

\n,l3;F,I,K,mF) = 

X] {j,mj,I,Tni\F,mF)\n;K,m,j)\/3;I,mi) , 

mj . rrij 

. (1) 

where n, k and j are the (Bohr) principal, the Dirac and 
the angular momentum quantum number of the electron 
respectively, while / represents the nuclear spin. Finally, 
mi, rrij and mp are the projections of the nuclear, elec- 
tronic and total (nucleus plus electron) angular momen- 
tum onto the quantization axis, respectively. Using stan- 
dard notation, {j,mj, I,mi\F,mp) are Clcbsch-Gordan 
coefficients. 

In writing Eq. ([T]) , we neglected any hyperfine interaction 
between nucleus and electron which could, in principle, 
be responsible of mixing ionic states with different quan- 
tum numbers. Such a mixing could in turn affect the 
angular and polarization properties of the emitted radi- 
ation, which are the physical quantities we are going to 
investigate. Since Eq. ([T|) will be used in what follows, 
we need to estimate this approximation. We may evalu- 
ate the goodness of Eq. ([1]) by calculating the admixture 
of the state nSi/2{F') that the state mSi/2{F) acquires, 
using the non-relativistic Hamiltonian which accounts for 
nucleus-electron hyperfine interaction |34| : 
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(2) 
where ag is the Bohr radius, kg = 2(77vMoA'-BMAf/47rft^ 
and Ho, Hb, Mjv, gN are the magnetic permeability of 
vacuum, the Bohr magneton, the nuclear magneton and 
the nuclear g-factor, respectively. In the last step of the 
above equation, the non-relativistic values of the ionic 
wavefunctions at the origin together with the operator 
relation 



Si= ifF2_j2_g2 



(3) 



have been used. Equation (jH]) is valid as long as the elec- 
tron states we use are spherically symmetric, that is if 
they are s-states. 



A property of Eq. ([2]) to remark is that only states with 
the same quantum number F can be coupled by the hy- 
perfine interaction Hamiltonian. For n = 1 and m ~ 2, 
the mixing coefficient C in Eq. ^ turns out to be of or- 
der 10^^ in Hydrogen atom while 10^^ in hydrogenlike 
Uranium ion. These two numbers can be taken as esti- 
mation for the accuracy of equation ([T]) , for the problem 
under consideration. 

Finally, due to parity invariance of the hyperfine interac- 
tion Hamiltonian, we underline that states with different 
parities may not get mixed. We will discuss our results in 
Sect, uni on the base of the considerations here presented. 
To further proceed, we suppose that the nucleus does not 
interact with the radiation field. In the language of quan- 
tum mechanics, this equates to considering that the inter- 
action Hamiltonian couples only electron fields through 
photon emission, while it does not act on the quantum 
space of nuclear states. This hypothesis holds for decays 
which involve ionic bound states, since the energy re- 
leased in such decays is far lower than the normal nuclear 
excitation energies. In contrast, in cases where the initial 
ionic state is a (positive) continuum-state, such hypoth- 
esis might fall as the energy of the radiation emitted in 
consequence of the electron capture might be enough to 
excite the nucleus and, therefore, to give rise to the so- 
called nuclear excitation by electron capture (NEEC) J27| . 
As a result of this supposition, we will find in the next 
subsection that the radial part of the decay amplitude 
will be characterized by only electron state components. 
Whereas, in the angular part, both the electron and nu- 
cleus states components will play a role, due to the cou- 
pling of the angular momenta. 



B. Second order transition amplitude 

The theory of two-photon decay is based on the sec- 
ond order transition amplitude, discussed for example in 
Akhiezer and Berestetskii [35|- For an initial state \i) and 
a final state |/), such amplitude reads 



M^^'^^i^ f) 
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where kj , wa ■ -^j and Uj are respectively the wave vec- 
tor, the polarization vector, the helicity and the energy 
of the jth emitted photon (j = 1,2). Ei^,^ arc the ener- 
gies of the initial and intermediate ionic state, while a 
is the standard vector of Dirac matrices. The subscripts 
i, I/, / will be used throughout the paper as referring 
to initial, intermediate and final state respectively. The 
summation over the intermediate states showed in for- 
mula ^ runs over the whole ionic spectrum, including 
a summation over the discrete part as well as an inte- 
gration over the (positive and negative) continuum. For 



the problem under consideration, sueh summation splits 
up into summations over the principal quantum number 
n^, the Dirac quantum number k^, the total angular mo- 
mentum F^ and its projection onto the quantization axis 

Since the two emitted photons define generally two differ- 
ent trajectories, the evaluation of the angular properties 
of the amplitude in Eq. Q is best carried out by de- 
composing the photon fields into their spherical tensor 
components of defined angular momentum. For the jth 
photon, such a decomposition reads [3a | 



WA,e' 



rtm) 



L = IM=-L (5) 



where [L] = 2L + 1 while -D^^ ^ . {kj — > z) are the Wigner 
rotation matrices which rotate each tensor component 



with angular momentum L and original quantization axis 
kj into the same component with the chosen quantiza- 
tion axis z. Furthermore, the standard notation A^^j^j 

and A^™/ is used to denote the electric and magnetic 
multipole, respectively. Each one of these multipoles can 
be expressed in terms of the spherical Bessel functions 
jhikr) and the spherical tensor Tlam as 
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Upon introducing Eq. ([S]) into (2]), expanding the ionic 
states as showed in Eq. ([T]) and by taking into account 
that the nuclear states must be normalized, the ampli- 
tude can be written as 
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where Aj runs from Lj — 1 to Lj -t- 1, T = i, 5 is used 
to denote the large (L) and small (S) components of the 
electron Dirac spinor, for which the factor P^ is defined 
as P^ — 1 and P^ = — 1, where T refers to the conjugate 
of T, i.e. T — L for T = S and vice versa. Furthermore, 
Pi^2 = 0, 1 and the function ^£^ is given by 
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for A = L - 1 

for A = L + 1 
otherwise . 



(8) 



The radial part of the amplitude in Eq. ([7]) is represented 
^"^"^A , which reads 

I Ao ' 



by the integral C/J"^ 



.,T' 



t^A/A2 =y drdr'r2r'2jAi(/cir')jA2(fer)g;*.gi;/+„^g: 

(9) 
where gj ^ are the small and large radial components of 
the final and initial electron state, while 



T T' 
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(10) 



is the radial Green function of the process. 

By inspecting Eq. ©, we notice that the radial part of 



the transition amplitude involves only electron state com- 
ponents. As mentioned in Sect. Ill Al this fact comes di- 
rectly from having neglected any photon-nucleus inter- 
action. Although the integral in Eq. ^ involves only 
electron state components, its evaluation is anyway a 
challenging task due to the (infinite) summation over 
the principal quantum number ny contained in the radial 
Green function. In the present work, such integral has 
been computed by using the Greens library 37| . This ap- 
proach has been presented for the first time in the context 
of two-photon decays by Surzhykov et al. in Ref. |14l | . 
The angular part of the amplitude in Eq. ([7]) is repre- 
sented by the elements xlai m and Xmi m ■ therein con- 
tained and can be computed analytically. With the help 
of Eq. ([Ij , we can indeed write 



x4, m,„ = E ^•?/ ' ™J/ J,mi\Ff,mFj) 



(11) 



1 fTT-F 



Xm, m,„ "= E ^^^ ' "*J» J,mi\F,,'i 

where a are Pauli matrices while the elements 
{Kf,lJ,mjf\a-fl^^^j^ijKi„l^,mjJ and {k„,1]^' ,mjja ■ 




in hydrogcnlike ions. It contains the complete informa- 
tion on the emitted radiation. For instance, the energy 
and angle-polarization distributions of the two emitted 
photons can be written in terms of such transition ampli- 
tude. Assuming that the ion is initially unpolarized and 
that the polarization of the final ionic state remain un- 
observed, the differential decay rate reads (atomic units) 



i 



dw 



Ai Ao 



UJ1UJ2 



dwidr^idr^a (27r)3c2 2F, + 1 



(13) 



FIG. 1: (Color online) Geometry considered for the two- 
photon emission. The propagation direction of the first pho- 
ton is adopted as z direction, x is chosen such that xz is the 
reaction plane (plane spanned by the two photon directions). 
6 is the angle between the photons directions, while angles 
Xi,2 define the linear polarizations of the first and second pho- 
ton respectively with respect to their respective polarization 
planes. The polarization plane of the first (second) photon is 
denoted by A (B) and represents the plane orthogonal to the 
photon direction. 
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iiAi 

where 



tdil 



If ,mj.) have been already discussed else- 
381 and will not be here recalled. 
The initial and final states involved in the two-photon 
transitions which we shall analyze in Sect. IIIII are 
unpolarized. Therefore, since for the decay of unpolar- 
ized states there is not any preferred direction for the 
overall system, we arbitrarily adopt the quantization axis 
(z) along the momentum of the "first" photon: z \\ ki. 
We furthermore adopt x such that the xz-plane is the re- 
action plane (plane spanned by the photons directions). 
Figure [T] sketches the geometry we consider for the decay. 
Within this geometry, the Wigner matrices in Eq. ([7]) 
simplify to 
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where d^ ^{6) is the reduced Wigner matrix and 9 is 
the polar angle of the second photon, which coincides, in 
the chosen geometry, with the angle between the photons 
directions (opening angle). Hence, the relative photons 
directions are uniquely determined by assigning the open- 
ing angle 9, which will be the variable against which the 
angular correlation will be plotted in Sect. IIIII 



C. Angular and Polarization correlations 



Within the approximations mentioned in Sect. Ill Al 
equation ^ represents the relativistic transition ampli- 
tude for the two-photon decay between hyperfine states 



By taking into account the axes geometry chosen for 
the two-photon emission, which has been explained in 
Sect. IIIB[ we can easily perform the integration over the 
solid angle of the first photon (di^i) as well as the in- 
tegration over the azimuthal angle of the second photon 
{d4>2)- Now, since part of this work is devoted to analyze 
photons polarizations, further details concerning the de- 
tection geometry must be provided, in order to proceed 
with the analysis 

In Fig. [1] we show how the photon polarizations may be 
defined in a case experiment. The polarization of each 
photon is to be measured in the "polarization plane" , 
which is the plane orthogonal to the photon direction. 
This ineans that, without loss of generality, we are con- 
sidering linear (not circular) polarizations of photons. 
This choice is motivated by the fact that, in high energy 
regime, linear polarizations of photons are much more 
easily measured than circular polarizations. In Fig. [TJ 
the polarization plane of the first and second photon are 
denoted by A and B, respectively. Each detector is sup- 
posed to have a transmission axis, along which the linear 
polarization of the photon is measured. Such a trans- 
mission axis is rotated by an angle x with respect to the 
reaction plane. Finally, each detector is supposed to work 
as a filter: Whenever a photon hits it, the detector gives 
off or not a "click" , which would respectively indicate 
that the photon has been measured as having its linear 
polarization along the direction x or x -f- 90°. 
By integrating over the first photon energy (dwi ) and by 
using the well-known relations between linear and circu- 
lar polarization bases [3^] 
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we can write the differential decay rate which depends 
upon the photons linear polarizations as 

^ ' dcos6i 2Fi + l ^ ^ 
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FIG. 2: (Color online) Angular correlation W{d) for the tran- 
sitions 2si/2 (F = 1,0) — >■ lsi/2 (F — 1,0) in Hydrogen atom. 
Tlie dashed-red curve refers to the electric dipole approxima- 
tion while the solid-black curve refers to the full multipole 
contribution. 



We shall call this function "polarization correlation". 
The polarization correlation represents the probability of 
detecting the emitted photons with defined linear polar- 
izations xi and X2, in the two-photon decays of hyperfine 
states in hydrogenlike ions. 

Finally, by summing over the photons polarizations, we 
define the "angular correlation" function as 



W{9) 



Stt^ 



dcos6l 2Fi + l 

mp^ mFfXiX2 



dcji 
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(16) 



Equipped with equations (|T5|) and ([T6| . in the next sec- 
tion we shall investigate the functions W^'^ ^'^ {9) and 
W{6) as obtained for a few decays of unpolarizcd hy- 
perfine 2s states in Hydrogen atom and in hydrogenlikc 
Uranium ion. 



III. RESULTS AND DISCUSSION 

In this section we shall analyze the angular and po- 
larization correlations defined in Eqs. P^ and p^ . for 
decays of hyperfine 2s states in Hydrogen atom as well 
as in hydrogenlike Uranium ion. 



A. Angular and polarization correlations in 
Hydrogen atom 

We first analyze Hydrogen atom, whose nuclear spin 
is, as known, I = 1/2. In particular, we analyze the four 
two-photon transitions 2si/2{F = 1,0) — )■ lsi/2 (-F — 



1,0). The predicted angular correlation W{9) for these 
transitions is displayed in Fig. [51 where the full mul- 
tipole and the electric dipole (ElEl) contributions are 
separately displayed. We might immediately notice that, 
in the two transitions 2si/2 (F = 1) — >■ lsi/2 {F = 1) 
and 2si/2 (-?" = 0) ^ lsi/2 {F = 0), the curves obtained 
within the electric dipole approximation practically co- 
incide with the curves obtained with the full multipole 
contribution. The angular correlation in these transi- 
tions can be well described by W{6) ~ 1 -I- cos^ 9, which 
is symmetric with respect to 9 = 90° and which is the 
analytical form that has been for long associated to the 
angular correlation in 2s — ^ Is two-photon transitions in 
Hydrogenlike bound systems [ll|, [13, |33] ■ 
Whereas, the radiation patterns in 2si/2 (F = 1) — >■ 
lsi/2 {F = 0) and 2si/2 {F ^ 0) ^ lsi/2 {F = 1) tran- 
sitions show a different, asym,m,etric shape. It can be 
seen that high-multipolc contributions are directly re- 
sponsible for the asymmetry, while the contribution of 
the leading ElEl multipole can be well described by 
Weiei{9) ^ 1 — 1/3 cos^ 9. This result is typical for two- 
photon decays of the type Jtot = 1 (0) — >■ Jtot = (1), 
where Jtot is the total angular momentum of the sys- 
tem which undergoes the decay. In fact the two-photon 
decay (ls2s)^S'j=i -^ (lsls)^S'j=o in heliumlike ions, 
where J is the total angular momentum of the system 
given by the coupling of the two electron angular mo- 
menta, shows approximately the same behavior jlTl |40| . 
Quantitatively, the asymmetry 
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1^(61 ==180°)- 14^(6* = 0°) 
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(17) 



is approximately 0.20 in F = 1 (0) ^^ F = (1) tran- 
sitions, while it is vanishing in F = 1 (0) ^ F = 1 (0) 
transitions. 

Difficulties in measuring the correlation function in those 
transitions where the total spin {F) is flipped (F = (1) 
— >■ F = 1 (0)) might arise from the fact that the de- 
cay rates for these transitions are noticeably nine orders 
of magnitude smaller than the decay rates for the tran- 
sitions in which the spin is not flipped (F = (1) — ?> 
F — (1)), as seen from Fig. [21 This fact is caused 
by the strong cancellation of the contributions given by 
the Pi/2 and P3/2 intermediate states (in the calcula- 
tion of the leading multipole ElEl) to the correlation 
function in F = (1) — > F = 1 (0) transitions. As a 
matter of fact, within non-relativistic dipole approxima- 
tion, these two contributions are equal and opposite, so 
that the transitions 2si/2 (-^ = 1) ~^ lsi/2 {F = 0) and 
2si/2 (^ = 0) ~^ 1*1/2 {F = 1) cannot proceed via two- 
photon dipole ElEl emission [I3, |4l|. Thus, the non- 
vanishing correlation functions showed in the top-right 
and bottom-left panels of Fig. [2] are only given by higli- 
multipole and relativistic effects. 

As a consequence of what said above, any experimental 
incoherent population of the initial and final hyperfine 
states would result in measuring the angular correlation 
function as given by ^ 1 -|- cos^ 9. This has been con- 
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FIG. 3: (Color online) Polarization correlation W^^ ^'■' (9) for 
the transitions 2si/2 (-F = 0) — !■ lsi/2 {F — 0) (solid-black 
curve) and 2si/2 {F — 0) -^ lsi/2 {F = 1) (dashed-green 
curve) in Hydrogen atom. The four polarization configura- 
tions (xi , X2 = 0° , 90° ) are separately displayed. 



firmed, for example, in Refs. [42|, |43[. For a direct ex- 
perimental investigation of the angular correlation func- 
tion in F = 1 (0) — )• F = (1) transitions, one would 
therefore normally need to selectively populate the ini- 
tial and observe the final hypcrfinc state. Such a selection 
could be made, in principle, by resolving energetically the 
emitted photons. The two hyperfine states 2s {F = 1) 
and 2s (F = 0), as well as the states Is (F = 1) and 
Is (F = 0), although they have been considered degener- 
ate in the present paper, possess in fact slightly different 
energies due to the well-known hyperfine energy split- 
ting, which is of the order ~ 10~^ eV in both cases (see 
Sect. Ill A|) . However, the energy resolution of normal pho- 
ton detectors for light quanta in the visible range, which 
is the energy range of the photons emitted in the con- 
sidered transitions, is far lower than 10~^ eV. Therefore, 
the selection of initial and final hyperfine states cannot 
proceed by analyzing energy. 

Alternatively, in order to achieve the same goal, we could 
make use of the linear polarization properties of the emit- 
ted photons. In Fig. [Sj we plot the polarization correla- 
tion W^'^ ^^{6) as obtained for the transitions 2si/2(F = 
0) ^ lsi/2(F = 0) and 2si/2(F = 0) ^ lsi/2(F = 1) 
in Hydrogen atom. As it can be seen, the orthogonal 
linear polarization configurations xi = 0°,X2 = 90° and 
Xi = 90°, X2 = 0° are absent inF=:O^F = transi- 
tion. More quantitatively, from Fig.|3]we can identify the 
photons' polarization state for F = — )■ F = transition 
as being 
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^=^==( 1X1=0, X2 = 0) 

\/l-cos^6'V (18) 

+ e'^cos0|xi = 9O°,X2 = 9O°)) , 



where the phase i5 can be determined, for instance, by 
analyzing the polarization correlation in the helicity ba- 
sis. By doing that, one easily finds (5 = 0. Although not 
directly shown, the photons polarization state in Eq. (J18p 
holds also in case the photons came from F = 1 — > F = 
1 transition. 

In Refs. [isllJli the same photons polarization state has 
been derived by analyzing 2s — > Is transitions in hydro- 
genlike ions within non-relativistic dipole approximation 
and without considering the nuclear spin. This is con- 
sistent with what found above, since, as mentioned, F 
= 0(1)— >F=1(0) transitions are forbidden within 
non-relativistic dipole approximation and therefore the 
polarization state displayed in Eq. (fT8|) is the only one 
which the photons can assume in two-photon 2s — > Is 
transitions. 

Coming back to Fig. [3l we notice that the decay rate 
for orthogonal polarization configurations is not vanish- 
ing for F = — >■ F = 1 transition, for any angle 6, 
although its order of magnitude is very small in com- 
parison with the decay rates shown in other panels of 
the same figure. This fact allows us to suggest a work- 
able scenario for investigating the angular and polariza- 
tion correlations in two-photon F = — > F = 1 transi- 
tions. By using standard devices of modern spectroscopy, 
we can experimentall y e xclusively populate the initial 
2si/2 (F = 0) state [4J]. Such a state would conse- 
quently decay either into lsi/2 (F = 0) or lsi/2 (F = 1) 
state. Now, as already underlined, only the transition 
2si/2 (^ = 0) ~^ lsi/2 {F = 1) of the two is character- 
ized by non- vanishing decay rate for two photons having 
orthogonal linear polarizations. Thus, we could select 
out the two-photon decay channel F = — > F = 1 just 
by filtering the photons polarizations. By virtue of the 
fact that the wavelength of each emitted quanta in this 
decay is roughly « 250 nm and that polarizer crystals 
arc available in such wavelength range, we could place 
one linear polarizer filter in front of each photon detec- 
tor. We would use the detectors in coincidence and we 
would set opposite linear polarizers transmission axes of 
the filters, in order to allow only the photon pairs com- 
ing from 2si/2 (^ = 0) — > lsi/2 (F = 1) transition to 
be detected. By inspecting the orders of magnitude of 
the decay rates in Fig. [31 we notice that, in order to suc- 
cessfully block the unwanted light from F = ^- F = 
transition, the extinction ratio of the polarizers must be 
higher than 10^'^:!. Such performances should be achiev- 
able by using in series, for instance, two or more normal 
polarizers whose extinction ratio is usually of order 10^:1. 
Problems related to the alignment of the polarizers might 
however arise. 

By using this sketched set-up, wc should therefore be sen- 
sible only to those photons that come from 2si/2 (F = 
0) — > lsi/2 (F = 1) transition and that have the se- 
lected polarization configuration, i.e. we will measure the 
dashed-green curve showed either in the bottom-left or in 
the top-right panel of Fig. [31 depending on the polarizers 
settings. 
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FIG. 4: (Color online) Angular correlation W(ff) for the tran- 
sitions 2si/2 (F = 4, 3) — ^ lsi/2 {F ~ 4, 3) in hydrogenlike 
92^U. The dashed-red curve refers to the electric dipole ap- 
proximation while the solid-black curve refers to the full mul- 
tipole contribution. 



Due to the intrinsic weakness of the signal coining from 
the desired decay channel, one could at this point think 
that the hyperfine interaction itself, here neglected from 
the outset, may mix states with different quantum num- 
bers, so that the angular and polarization properties of 
the emitted radiation could vary from what shown in 
Figs. [2] and [H However, as briefly discussed in Sect. Ill A"l 
the Hamiltonian which accounts for the hyperfine inter- 
action neither mixes states with different parities nor S*- 
states with different quantum number F. This implies 
that the state 2si/2 (F = 0) docs not acquire any admix- 
ture of any other state and thus remains well defined. 
Indeed low energy hyperfine states in Hydrogen atom are 
known to be well separated, as it is evident from the fact 
that the light emitted in the hyperfine transition of Isi /2 
state in Hydrogen is one of the most accurately measured 
physical quantities ever and is of great importance in as- 
trophysics and cosmology. 



B. Angular and polarization correlations in 
hydrogenlike Uranium 

We proceed to analyze hydrogenlike g|^U ion, whose 
nuclear spin is I = 7/2 |45|. In particular, we study the 
transitions 2si/2 {F = 4, 3) ^ lsi/2 {F = 4, 3) in such 
ion. Hydrogenlike heavy ions are interesting bound sys- 
tems both from a theoretical and an experimental point 
of view, since they are as simple as Hydrogen atom but 
show, at the same time, remarkable relativistic and high- 
multipole effects. We have nonetheless seen that high- 
multipolc effects show already upinF = 0(l)^-F = l 
(0) transitions in Hydrogen atom, where they are directly 
responsible for the asymmetric behaviour of the angular 



correlation function (cfr. Fig. [2]). In |2^U ion, we expect 
an enhancement of this effect. 

The angular correlation, W{6), of the photon pair for 
2si/2 {F = 4, 3) -> lsi/2 {F = 4, 3) two-photon tran- 
sitions in hydrogenlike 92^ U ion is displayed in Fig. |4l 
The full multipole and the electric dipole (ElEl) con- 
tributions are separately displayed. Although the spin 
properties of |2^U and JH are considerably different, we 
can see that the relative (photon) angular correlations 
arc rather similar in shape. However, by inspecting care- 
fully both Figs, m and [21 it can be noticed that, as ex- 
pected, high-multipolc effects play overall a more impor- 
tant role in g2^U than in Hydrogen. First, in contrast 
to Hydrogen, the angular correlation for non-spin-flip 
transitions (F = 3 (4) ^- F = 3 (4)) show sizable de- 
viations from the ~ 1 + cos^ shape. This effect is al- 
ready known from the past literature, where it has been 
showed that high multipoles contribute with terms of the 
type '^ cos 9 to the angular correlation in 2si/2 ^ lsi/2 
transitions [l9|. Secondly, we also notice that, in spin- 
flip transitions (F = 3 (4) — >■ F = 4 (3)), deviations of 
the angular correlation function from the symmetric an- 
alytical formula ~ 1 — l/3cos^0 are more pronounced 
in hydrogenlike g^^U ion than in Hydrogen. More con- 
cretely, for the case of 92^!! ion, the asymmetry A de- 
fined in Eq. ([T7]) turns out to be approximately 0.26 in 
2si/2 (-F = 4(3)) -^ lsi/2 {F = 3(4)) transitions, while 
it is 0.09 in 2si/2 {F = 4(3)) -^ lsi/2 {F ^ 4(3)) transi- 
tions. 

It can be furthermore noticed that each one of the 2s — ?> 
Is two-photon transitions shown in Fig. 3] is character- 
ized by far higher decay rate with respect to the transi- 
tions shown in Fig. [21 This fact is not surprising since 
it is known that the total decay rate for this transition 
rises fast with Z, as it can be approximately expressed 
by crw 8.226 Z^ sec'^ 0. 

As a final remark, by comparing Figs. [H and [H we no- 
tice that spin-flip transitions are less suppressed in hy- 
drogenlike Uranium than in Hydrogen atom, which is a 
consequence of the fact that relativistic and retardation 
effects characterize highly charged ions. The selective de- 
tection of the photon pair coming from F = 3 (4) ^ F 
= 4 (3) transitions in hydrogcnlike Uranium, however, is 
complicated by the fact that the energy of the emitted 
radiation is in the range of hard X-rays. The discernment 
by energy analysis of the decay channels in g2^U share the 
same problem already encountered in the Hydrogen case: 
The hyperfine energy splitting between 2s]^ ^ {F — 4) and 
2si/2 {F — 3), as well as the one between lsi/2 {F = 4) 
and lsi/2 {F ~ 3), is of the order ^ 1 cV, which is far 
smaller than the normal resolution of X-ray detectors. In 
addition, since there exist no polarizer filters for such en- 
ergy regime, the photon polarization analysis suggested 
to select the desired decay channel in the Hydrogen case 
is not immediately applicable to Uranium, regardless of 
the polarization properties of the emitted photons. As 
a matter of fact, the polarization resolved experiments 
in X-ray energy regimes are not normally performed by 
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FIG. 5: (Color online) Polarization correlation W^^ ^^ (9) for 
the transitions 2si/2 (-F = 3) — ^ lsi/2 {F = 3) (solid-black 
curve) and 2si/2 {F — 3) ^ lsi/2 {F = 4) (dashed-green 
curve) in hydrogenlike 92^ U ion. The four polarization con- 
figurations (xiiXa = 0°,90°) are separately displayed. 



using polarizer filters, but rather by using Compton po- 
larimeters |28l - l3(]l l46l l47l|. In such devices, the informa- 
tion about the photon polarization is obtained from the 
reconstruction of the Compton scattering events (of the 
incident photons) in the detector and thus has a statisti- 
cal nature. As a consequence, Compton polarimeters are 
not normally used to record the polarization properties 
in multi-photon processes. Rather, they are employed to 
record the polarization properties in single photon pro- 
cesses [3l|. Nonetheless, by selecting events which have 
been recorded in coincidence by two or more Compton 
polarimeters and which have the desired scattering an- 
gle, information on the polarization of a multi-photon 
state could be in principle achieved [26, 48 51]. 
In the light of such possibility, we present in Fig. [5] 
the function W^'^^^{9) as obtained for the transitions 
2si/2(^ = 3) ^ I.si/2(F ^ 3) and 2si/2(F = 3) ^ 
lsi/2(-F = 4) in hydrogenlike g2^U ion. We see from the 
figure that the general scheme for which photons com- 
ing from spin-flip transitions (non-splin flip transitions) 
have mainly orthogonal (parallel) linear polarizations re- 
mains valid also in hydrogenlike Uranium. However, in 



the Uranium case, the suppressed polarization configura- 
tions are more evident, especially in spin-flip transitions, 
where they turn out to be of the same order of magni- 
tude of the leading polarization configurations. Thus, in 
contrast to Hydrogen, there exist not any polarization 
configuration by which to single out the photons coming 
from one transition of the two displayed in Fig. \E\ 
IV. SUMMARY 

In summary, the angular and polarization correlations 
in two-photon decays of 2s hyperfine states in hydro- 
genlike ions have been investigated within second-order 
perturbation theory and the Dirac relativistic framework. 
Results for the angular correlation have been showed for 
the transitions 2si/2 (F = 1, 0) ^ lsi/2 (F = 1, 0) in Hy- 
drogen atom as well as for 2si/2 {F = 4, 3) — > lsi/2 {F = 
4, 3) in hydrogenlike |2''U ion. It has been possible to 
identify two types of emission pattern: 1- the two-photon 
decays which connect ionic states with the same spin 
(non-spin-flip transitions) are found to be characterized 
by an angular correlation W{9) '--^ 1 -I- cos^ 9, with small 
deviations which are increasing with the atomic number 
Z; 2- the two-photon decays which connect ionic states 
with different spin (spin-flip transitions) are found to 
have much smaller decay rate and to be characterized by 
an angular correlation W{9) ^--^ 1 — 1/3 cos^ 6, with devia- 
tions which are already important for low-Z systems. In 
both cases, the (asymmetric) deviations come from high- 
multipole contributions. 

Results for the polarization correlation have been showed 
for the transitions 2si/2 (^ = 0) — >■ lsi/2 {F = 1,0) in 
Hydrogen as well as for 2si^2 (^ = 3) — ^ lsi/2 {F = 4, 3) 
in hydrogenlike g2^U ion. It has been found that the 
photons coming from the studied F = ^ F = tran- 
sition in Hydrogen atom are exclusively characterized by 
parallel linear polarization. By using such results, it has 
been possible to outline a way to experimentally measure 
the angular and polarization correlations in the transition 
2si/2 (F = 0) — > lsi/2 {F = 1), in Hydrogen atom. 
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